Proyecto clustering de imágenes

Autor: Josué Fabricio Urbina González

Se utilizaron los siguientes módulos de python para su realización.


In [1]:
%matplotlib inline
from tensorflow.examples.tutorials.mnist import input_data
from sklearn.manifold import TSNE
from sklearn.decomposition import PCA
from sklearn.cluster import AgglomerativeClustering
import numpy as np
import matplotlib.pyplot as plt

In [ ]:
sns.set_style('darkgrid')
sns.set_palette('muted')
sns.set_context("notebook", font_scale=1.5,
                rc={"lines.linewidth": 2.5})

Se leen los datos del Fashion-MNIST. Y se declara un RandomState utilizado para el TSNE.


In [ ]:
# Visual
# Random state.
RS = 102030

data = input_data.read_data_sets('Data/P2')

tmp = data.train.images
tmp_labels = data.train.labels

Reordenamos los puntos de acuerdo a las clases con las que fueron etiquetadas.


In [ ]:
X = np.vstack([tmp[tmp_labels == i]
               for i in range(10)])
y = np.hstack([tmp_labels[tmp_labels == i]
               for i in range(10)])

Seleccionamos un conjunto aleatorio de 6000 muestras.


In [ ]:
set = set()
while len(set) < 6000:
    set.add(np.random.randint(0, len(tmp)))
set = sorted(list(set))
train = X[list(set), :]
train_labels = y[list(set)]

Realizamos la reducción de dimensionalidad a 2 para su visualización. Utilizando TSNE


In [ ]:
proj = TSNE(random_state=RS, perplexity=50, verbose=1).fit_transform(train)

La siguiente función es de apoyo para mostrar los puntos


In [ ]:
def scatter(x, colors):
    palette = np.array(sns.color_palette("hls", 10))

    # We create a scatter plot.
    f = plt.figure(figsize=(8, 8))
    ax = plt.subplot(aspect='equal')
    sc = ax.scatter(x[:,0], x[:,1], lw=0, s=40,
                    c=palette[colors.astype(np.int)])

    plt.xlim(-25, 25)
    plt.ylim(-25, 25)
    ax.axis('off')
    ax.axis('tight')

    return f, ax, sc

Usando reducción mediante PCA y el método AgglomerativeClustering se tratan los datos.


In [ ]:
pca = PCA(n_components=300)
proj = pca.fit_transform(train)

# Agglomerative Clustering
cluster = AgglomerativeClustering(n_clusters=10)
cluster.fit(proj)

In [ ]:
scatter(proj_v, cluster.labels_)
plt.savefig('fashion_tsne-generated.png', dpi=120)
plt.show()

In [ ]:
scatter(proj_v, train_labels)
plt.savefig('fashion_tsne-generated_original.png', dpi=120)
plt.show()